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Abstract 



The presence of chaos and quantum chaos is shown in two different nuclear 
systems. We analyze the chaotic behaviour of the classical SU(2) Yang-Mills- 
Higgs system, and then we study quantum chaos in the nuclear shell model 
calculating the spectral statistics of A = 46-50 atomic nuclei. 



1. Introduction 

A classical system is chaotic if has exponential divergence of initially 
closed trajectories in the phase space, i.e. small differences in the initial con- 
ditions produce great changes in the final behaviour. The phenomenon of 
chaos, which has become very popular, rejuvenated interest in nonlinear dy- 
namics. Through numerical simulations on modern computers the presence 
of chaos has been discovered to be pervasive in many dynamical systems of 
physical interest [1]. 

In quantum mechanics one cannot apply classical concepts and methods 
directly being the notion of trajectory absent. Nevertheless, many efforts 
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have been made to establish the features of quantum systems which reflect 
the qualitative difference in the behaviour of their classical counterparts [2]. 

The quantum chaos, or better the quantum chaology [3] , studies the prop- 
erties of quantum systems which are classically chaotic. The energy fluctu- 
ation properties of systems with underlying classical chaotic behaviour and 
time-reversal symmetry agree with the predictions of the Gaussian Orthog- 
onal Ensemble (GOE) of random matrix theory, whereas quantum analogs 
of classically integrable systems display the characteristics of the Poisson 
statistics [2-4]. 

In this paper we analyze the presence of chaos and quantum chaos in two 
systems of interest for nuclear physics. In the first part we study analyt- 
ically the suppression of classical chaos in the spatially homogenous SU(2) 
Yang-Mills-Higgs (YMH) system. In the second part we analyze numerically 
quantum chaos in the nuclear shell model by studying the spectral statistics 
of low-lying states oi A — 46-50 atomic nuclei. 

2. Chaos in the SU(2) Yang— Mills— Higgs system 

In the last years there has been much interest in the chaotic behaviour of 
classical field theories [5] . Usually the order-chaos transition in these systems 
has been studied numerically by using Lyapunov exponents and Poincare sec- 
tions [1,2]. Less attention has been paid to analytical criteria [6,7]. Here we 
study the spatially homogenous SU(2) YMH system. Obviously, the constant 
field approximation implies that our SU(2) YMH system is a toy model for 
classical non-linear dynamics, with the attractive feature that the model 
emerges from particle physics. 

The lagrangian density for the SU(2) YMH system is given by [8]: 

L = -If^^^F'^^'^ + \{D,<t>nD>^<t>) - V{<t>) , (1) 

where: 

F;, = d,Al - d^Al + ge^^^A^Al , (2) 

{D,4>) = d,4>-tgAlT'(t>, (3) 

with = b = 1,2,3, generators of the SU(2) algebra, and where the 

potential of the scalar field (the Higgs field) is: 

v(0)=/.2|0p+Ai0r (4) 
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We work in the (2+l)-dimensional Minkowski space (/i = 0, 1, 2) and choose 
spatially homogenous Yang-Mills and the Higgs fields: 

^^^^ = 9^0 = 0, i = l,2 (5) 

i.e. we consider the system in the region in which space fluctuations of flelds 
are negligible compared to their time fluctuations. 

In the gauge — and using the real triplet representation for the Higgs 
field we obtain: 

L = + io + ^ - g'illlAl - • 1,)'+ 

+ (Al + Al)$' - {A ■ $f - (^2 • $)^\ - V{$1 (6) 

where 0= (0S 02,03), Ai = {A\,AlA^), and A2 = {AlAlA^). 

When //^ > the potential V has a minimum in |0| = 0, but for //^ < 
the minimum is: 

which is the non zero Higgs vacuum. This vacuum is degenerate and after 
spontaneous symmetry breaking the physical vacuum can be chosen 0o — 
(0,0,^;). If A\ — Qi, Al — q2 and the other components of the Yang-Mills 
fields are zero. The system hamiltonian in the Higgs vacuum is: 

H-l{pi+ pI) + g'v^ql + ql) + \g'qlql , (7) 

where pi = qi and p2 = q2- Obviously — 2g^v^ is the mass term of the 
Yang-MiUs fields. 

The transition order-chaos in systems with two degrees of freedom may 
be studied by the curvature criterion of potential energy [6]. It is however 
important to point out that in general the curvature criterion guarantees 
only a local instability and should therefore be combined with the Poincare 
sections [9]. 

At low energy the motion near the minimum of the potential: 

V{qi, q2) = g'v^ql + ql) + \g''qlql (8) 
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where the curvature is positive, is periodic or quasiperiodic and is separated 
from the instabihty region by a hne of zero curvature; if the energy is in- 
creased, the system will be for some initial conditions in a region of negative 
curvature, where the motion is chaotic. According to this scenario, the en- 
ergy of order^chaos transition is equal to the minimum value of the line 
of zero gaussian curvature K{qi,q2) on the potential-energy surface. For 
our potential the gaussian curvature vanishes at the points that satisfy the 
equation: 

m - = + '''^^^">''' + - ^ 

It is easy to show that the minimal energy on the zero-curvature line is given 
by: 

= Vrmn{K = 0, Qi) = 6g^v\ (10) 

and occurs at qi = ±\/2v. Therefore the chaos-order transition depends on 
the parameter e = g^v^/ E: for < e < 6 a relevant region of the phase-space 
is chaotic, while for e > 6 the system becomes regular. This result shows 
that the Higgs field value in the vacuum v plays an important role: for large 
values it makes the system regular in agreement with previous numerical 
calculations of Savvidy [5]. Also the Yang-Mills coupling constant g has 
the same role. Instead for fixed v and g there is an order-chaos transition 
increasing the energy E [10]. 



3. Quantum chaos in the nucleeir shell model 



One of the best systems for the study of quantum chaos is the atomic 
nucleus, which has been the subject of many investigations [11]. In atomic 
nuclei, the fluctuation properties of experimental energy levels are best stud- 
ied in the domain of neutron and proton resonances near the nuclcon emission 
threshold, where a large number of levels with the same spin and parity in the 
same nucleus are present, and an excellent agreement with GOE predictions 
has been found [12]. 

In this work we undertake the statistical analysis of the shell-model en- 
ergy levels in the A = 46-50 region. By using second-quantization notation, 
the nuclear shell-model hamiltonian may be written as [13]: 

H ^^eaa^aa+ < 0(f3\V\5j > a^a'^ajUs, (11) 
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where the labels denote the accessible single-particle states, is the corre- 
sponding single-particle energy, and < q;/9|V"|57 > is the two-body matrix 
element of the nuclear residual interaction. 

Exact calculations are performed with the /7/2, /5/2) andpi/2 single- 
particle states, assuming a ^°Ca inert core. The diagonalizations are per- 
formed in the m-scheme using a fast implementation of the Lanczos algo- 
rithm with the code ANTOINE [14] . For a fixed number of valence protons 
and neutrons we calculate the energy spectrum for projected total angular 
momentum J and total isospin T . The interaction we use is a minimally 
modified Kuo-Brown realistic force with monopole improvements [15]. 

We calculate the T = states from J = to J = 9 for all the combi- 
nations of 6 active nucleons, i.e. ^^V, ^^Ti, ^^Sc, and "^^Ca, and also for ''^Ca 
and ™Ca. 

Since we are looking for deviations from chaotic features, we are mainly 
interested in the low-lying levels up to a few MeV above the JT yrast line. 
For each JT set of levels the spectrum is mapped into unfolded levels with 
quasi-uniform level density by using the constant temperature formula [16]. 

The spectral statistic P{s) is used to study the local fluctuations of the 
energy levels. -P(s) is the distribution of nearest-neighbour spacings si — 
{Ei^i — Ei) of the unfolded levels Ei. It is obtained by accumulating the 
number of spacings that lie within the bin (s, s + As) and then normalizing 
P{s) to unity. 

For quantum systems whose classical analogs are integrable, P{s) is ex- 
pected to follow the Poisson limit, i.e. P{s) — exp (— s). On the other hand, 
quantal analogs of chaotic systems exhibit the spectral properties of GOE 
with P{s) = (7r/2)sexp(-fs2) [2-4]. 

The distribution P{s) is the best spectral statistic to analyze shorter 
series of energy levels and the intermediate regions between order and chaos. 
The P{s) distribution can be compared to the Brody distribution [17] 

P{s,uj) =a{(jj + l)s'^exp{-as'^+^), (12) 

with 

a^{T[^]r\ (13) 

This distribution interpolates between the Poisson distribution {lo — 0) of 
integrable systems and the GOE distribution {lo — 1) of chaotic ones, and 
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thus the parameter uj can be used as a simple quantitative measure of the 
degree of chaoticity. 

In order to obtain more meaningful statistics, P{s) is calculated using 
the unfolded level spacings of the whole set of J = 0-9 levels for fixed T up 
to a given energy limit above the yrast line. Thus the number of spacings 
included is reasonably large. 



Energy 


46V 




4«Sc 


^«Ca 


4«Ca 


5"Ca 


< 4 MeV 


1.14 


0.90 


0.81 


0.41 


0.58 


0.67 


< 5 MeV 


1.10 


0.81 


0.96 


0.53 


0.58 


0.69 


< 6 MeV 


0.93 


0.94 


0.99 


0.51 


0.66 


0.62 



Table 1: Brody parameter uj for the nearest neighbour level spacings distri- 
bution for < J < 9, T = states up to 4, 5 and 6 MeV above the yrast 
line in the analyzed nuclei. 

Table 1 shows the best fit Brody parameter oj of the Pis) distribution for 
the J = 0-9 set of level spacings in the A = 46 nuclei up to 4, 5 and 6 MeV 
above the yrast line. Clearly, ^^V, ^^Ti and ^^Sc arc chaotic for these low 
energy levels, but there is a considerable deviation from GOE predictions in 
^^Ca, which is a single closed-shell nucleus. In view of the peculiarity of this 
nucleus, we performed calculations for ^^Ca and ^'^Ca, and obtained again 
strong deviations toward regularity, as the values of Table 1 show. 

To explain these results we observe that the two-body matrix elements 
of the proton-neutron interaction are, on average, larger than those of the 
proton-proton and neutron-neutron interactions. Consequently the single- 
particle mean-field motion in nuclei with both protons and neutrons in the 
valence orbits suffers more disturbance and is thus more chaotic. 

4. Conclusions 

The presence of chaos in the homogenous SU(2) Yang-Mills-Higgs system 
has been studied. The results show that for large values of the Higgs field in 
the vacuum there is a suppression of the chaotic behaviour of the Yang-Mills 
fields. In the future it will be of great interest to analyze this effect also in a 
non homogenous situation. 
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In the study of the nuclear shell model we have seen that the disturbance 
of single-particle motion by the two-body interaction is greatest in light 
nuclei, where the size of the range of the single-particle orbits is not much 
longer than the range of the nuclear force. In particular, for Ca isotopes, we 
find significant deviations from the predictions of the random-matrix theory 
which suggest that some spherical nuclei are not as chaotic in nature as the 
conventional view assumes. 
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